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We present a new independent evaluation of the hadronic vacuum polarization contributions aJ lad (Z.o) = 
§(<? — 2)'l ad (l.o) to the anomalous magnetic moment {anomaly) of the muon and tau leptons using r-decays 
and e + e~ data. The alone theoretical input used for describing the high-energy region not accessible exper- 
imentally is perturbative QCD plus (negligible) additional effects due to the QCD vacuum condensates. We 

11 and a^ ad (l.o) = 353.6(3.8) x 1(T 8 , which we compare with previous 
determinations. Our analysis leads to the Standard Model (SM) prediction: a^ M = 116 591 846.9(78.9) x 1CT 11 . 



obtain: a^ ad (l.o) = 7020.6(75.6) x 10~ 



Confronting 



with the recent BNL measurement leads to a^f w 

had { 



176(170) x 10" 



Combined 
li 



with the mean existing determinations of a^ ad (Lo), this leads to the conservative range: -42 < a" em x 10 11 < 413 
at 90% CL, from which we derive lower bounds on the scales of some new physics. We also update our 
old and first predictions of the SM contributions to a T . Including QED to sixth order, higher order hadronic 
and electroweak contributions, we obtain af M — 117 755(7) x 10 -8 , waiting for a future precise measurement of a T . 
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1. INTRODUCTION 

The recent E821 BNL result for the measurement of 
the muon anomaly averaged with older determinations 
[11 gives §: 



exp 



116 592 023(151) x 10 



-li 



(1) 



The announcement of a 2.6a deviation |lJJ4] of this ex- 
perimental result from the standard model (SM) pre- 
diction has stimulated intensive publications on the 
possible signal or/and constraints on new physics, but 
has also raised some criticisms ||. In this paper, 
we consider that the only rational answer to the the 
criticisms of [Bj is an independent re-evaluation of 
the hadronic vacuum polarisation contribution to the 
muon anomaly as an update of our old estimate in , 
which is the main motivation of this work. In order 
to avoid some specific theoretical dependences of the 
result, we shall mainly use the available data from t- 
decay and e + e~ — > hadrons, and limit ourselves to the 
use of perturbative QCD plus (negligible) additional 



effects due to the QCD vacuum condensates for de- 
scribing the high-energy region (QCD continuum) not 
accessible experimentally. We shall compare our re- 
sults with previous determinations and analyze its im- 
plications for some new physics (e.g. supersymmctry, 
radiative muon, composite, extended Zee and lepto- 
quark models) beyond the SM. Finally, we update our 
old and first |Q SM predictions for a T . 

2. THE HADRONIC VACUUM POLARISA- 
TION CONTRIBUTION TO THE LEP- 
TON ANOMALY 

Several papers have been devoted to the analysis of 
the contribution to the muon anomaly. Works prior 76 
have been reviewed in jq , while more recent works (af- 
ter 85) are reviewed in A partial historical review 
of the different determinations since 61 can be found 
in Figure 2 of ||, which we complete and update in 
Table 1, where only works published in journals with 
referee-policies have been considered. Using a disper- 
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Table 1 

Time and precision evolutions of the determinations of a^ ad (/.o) from e + e~ data. Results marked with * are our 
selected results which take into account update and results with minimal theoretical inputs. 



a£ ad (l.o) x 10 11 Authors Ref. Comments 



« 3200 

5500 - 11000 

fa 7500 

34 00+ 1900 

6500(500) 

6800(900) 

7300(1000) 

6630(850) 

6990(880) 

7020(800)* 

6840(110) 

7068(174)* 

7100(116)* 

7050(78) 

7024(153)* 

7113(103) 

7026(160)* 

6950(150)* 

7011(94) 

6951(75)* 

6924(62) 

7021(76)* 



Bouchiat-Michel (61) 
Durand (62) 
Kinoshita-Oakes (67) 
Bowcock (68) 
Gourdin-de Rafael (69) 
Bramon-Etim-Greco (72) 
Bailey et al. (75) 
Barger-Long-Olsson (75) 
Calmet-Narison-Perrottet-de Rafael (76&77) 
Narison (78) 
Barkov et al. (85) 
Kinoshita-Nizic-Okamoto (85) 
Casa-Lopez-Yndurain (85) 
Dubnicka-Martinovic (90) 
Eidelman-Jegerlehner (95) 
Adel-Yndurain (95) 
Brown- Worstell (96) 
Alemany-Davier-Hocker (97) 
Alemany-Davier-Hocker (97) 
Davier-Hocker (98) 
Davier-Hocker (98) 

Narison (01) 



ft 

[10] 

1 11 1 

[12! 

[13! 

\u\ 

[15[ 
[16[ 

if 

|7| 

[17[ 

[18[ 

[19! 

[20] 

[21! 
[22] 

[23] 
[24] 
[24] 
[25] 
[26] 



This work 



7T ' 7T 
7T + 7r" 



only 
only 



+r-decay 
+r-decay 

+Thcorctical input 

Average from e + e~ and 
r-decay (see section 3) 



Informative weighted mean values 



6992(28) stat 
7008(41) rtat 



Since 69 

Selected results marked with * 



sion relation, the hadronic vacuum polarisation contri- 
bution to the lepton anomaly ai can be expressed as 



1 Z" 00 

a? ad (l.o) = — dt K,(t) a H (t) . 
— Ki{t> 0) is the QED kernel function 



1 , x 2 {l-x) 
ax- 



with the analytic forms: 

Xi(t<4m?) = ~-4 yi - 

4yi{l-2yi)log(4yi) 
-2 (1 - 8yi + 8yf) x 

Vvii 1 - Vl) arccos y/yi 



(2) 



(3) 



zi 



with: 



yi 



4mf ' 



zi 



(1 + vi) 



l-zi 



and vi 



log(l + zi) 

2? log zi, 



1 



4mf 
t 



(4) 



(5) 



Ki{t) is a monotonically decreasing function of t. For 
large t, it behaves as: 



(G) 



which will be useful for the analysis in the large t 
regime. Such properties then emphasize the impor- 
tance of the low-energy contribution to a^ ad (l.o) (I = 
e, /i), where the QCD analytic calculations cannot be 
applied. 

— <?\tr(i) = c(e + e~ — * hadrons) is the e + e~ — > hadrons 
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total cross-section which can be related to the hadronic 
two-point spectral function Imtl(i) em through the op- 
tical theorem: 

„ a(e + e~ — > hadrons) 

R e +e- = -V+ t r = 127rlmn (t ) em , 7 

where: 

<r(e + e~ — > p + p~) 
Here, 



is the charged vector two-point correlator: 



47TCT 

3* 



(8) 



IC; = ijd^x e«* (0\TJ^(x) (J^(i)) t |0) 

= -(.^Y-^)iW<z 2 ) (9) 

is the correlator built from the local electromagnetic 
current: 

= - |dy*d - 5*7"* + ... (io) 

— In the following, we shall discuss in details the dif- 
ferent hadronic contributions to a t Jf d {l.o). Our results 
will be given in Table 2. 

3. HADRONIC VACUUM POLARISA- 
TION CONTRIBUTION TO THE MUON 
ANOMALY 

3.1. LIGHT ISOVECTOR MESONS 
Region below 0.8 GeV 2 

Due to the i-behaviour of the kernel function K^{t), 
this region gives so far the most important contribu- 
tion to a^ ad (l.o) (w 68% of the total contribution), and 
also the largest source of the errors (77% when added 
quadratically) . In our numerical analysis, we divide 
this region into three subregions. 

— The first one is the region near the two pions thresh- 
old 4m 2 < t < 0.4 GeV 2 , where the pion form factor 
l-FVl 2 ^) is constrained by universality |F T | 2 (0) = 1 and 
by the t expansion predicted by chiral perturbation 
theory: 



\F^\ 2 (t)^l + ^(r 2 )^ + C7r t 2 + 0(t 3 ) , 



(11) 
: (0.431 ± 

0.026) fm 2 is the mean pion charge radius squared 
and a* ~ (3.2 ± 1.0) GeV 4 pi. One can inspect that 



or its resummed expression 



(r 2 



the ALEPH/OPAL gfU) and e+e" data compiled in 
pl^pip^ ] satisfy both constraints. 

— The second region is the one from 0.4 GeV 2 < t < 
0.6 GeV 2 on the top of the p resonance. 

— For these two regions, we shall use CVC hypoth- 
esis which relates the electromagnetic to the charged 
current through an isospin rotation p2|] : 



. . 4-7ra 2 

(TH{t) = —j— v l ■ 

We follow the notation of ALEPH 



where: 



ImiP 



1 

2^r' 



(12) 



(13) 



ll ud,V 



i / fx e ^(o|r^(x)(j^(o)) T |o) 



(14) 



built from the local charged current J£dv( x ) = 
uj^dfa) . We use the accurate semi-inclusive ALEPH/ 
OPAL data [ j30|]3"l|l . For a comparison, we also show 



the results obtained from the use of ■ 



data in the 



whole region below 0.6 GeV 2 . 

— The third region is the one from 0.6 GeV 2 to 0.8 
GeV 2 , which is peculiar due to the uj — p mixing. In this 
region, we use either the r-decay data with corrected 
uj — p mixing effect, or the e + e~ data alone compiled 
in |p4 30 1 . One can check from the data that the effect 
of the 47T is negligible in such a region. 

Region from 0.8 to 3 GeV 2 

For our fitting procedure, we use CVC and the 
ALEPH/OPAL semi-inclusive r-decay data [|o|j3l). 
The dominant error from this region comes from the 
one between 2.1 and 3 GeV 2 (88% when the errors are 
added quadratically), due to the incaccuracy of the 
data near M T . For a comparison, we also show the 
results when using the sum of exclusive modes from 
e + e~ compiled in ]30p^| , where in this case the errors 
come mainly from the region below 2 GeV 2 . 

3.2. LIGHT ISOSCALAR MESONS 

u> and mesons 

We treat these mesons in a narrow width approxi- 
mation (NWA), which is expected to be a good ap- 
proximation. Using the relation between the e + e~ — > 
hadrons total cross section and the leptonic width T ee 
of the meson with a mass Mr: 



a H {t) NWA ~ 127r 2 ^g<5(t - M 2 ) 



one can write: 



had 



{1.0) 



3 r e 



NWA 



(15) 



(16) 



7tM R 

where, we shall use the PDG values of the electronic 
widths H). 

Region below 1.39 2 GeV 2 

We estimate the effect of this region by using the sum 
of exclusive 1 = modes for the ratio °_ compiled 
in@. 

Region between 1.39 2 and 3 GeV 2 
In order to account for the multi-odd pions, 
K K, KKir, ... modes, we estimate the effect of this 
region by assuming that it is mainly given by the 
w(1419), w(1662) and 0(1680), with their parame- 
ters measured by the DM2 Q and DM1 Q col- 
laborations. We shall estimate their leptonic widths 
which are multiplied by the hadronic branching ratios 
in ||^,|3(|, by assuming (see PDG and the indi- 
cation from DM2 and DM1 data), that the w(1419) 
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Table 2 

Determinations of a^ ad (l.o) using combined e + e~ and inclusive r decay data (2nd and 4th columns) and averaged 
e + e~ data (3rd column). 



Region in GeV 2 



a£ ad (l.o) x 10 1 



a? ad (l.o) x 10 s 



Data input 



Light Isovector 

Ami 0.8 
0.8 -> 2.1 
2.1 3. 
Total Light 1=1 
Light Isoscalar 
Below 1.93 



r+e + e 

4794.6 ±60.7 
494.9 ± 15.8 
202.0 ±29.7 
5491.5 ± 694 



4730.2 ±99.9 
565.0 ±54.0 
175.9 ± 16.0 
5471.1 ± 114. 



r+e + e 

165.8 ± 1.5 
28. 7± 1.1 
17.0 ±2.6 
211.5 ± 3.2 





387.5 ± 13 


387.5 ± 13 


15.3 ±0.5 





393.3 ±9.9 


393.3 ±9.9 


21.0 ±0.5 


0.66 -> 1.93 


79.3 ± 18.7 


79.3 ± 18.7 


4.3 ± 1.1 


From 1.93 to 3 








w(1.42), w(1.65) 


31.3 ±6.8 


31.3 ±6.8 


2.6 ±0.7 


0(1.68) 


42.4 ± 18.2 


42.4 ± 18.2 


3.8 ± 1.3 


Toted Light L=0 


933.5 ± 31.5 


933.S ± 31.5 


47.0 ± 2.0 


Heavy Isoscalar 








J/V(15 -» 4.415) 


87.0 ±4.7 


87.0 ±4.7 


13.08 ±0.69 


Y(15^ 11.020) 


0.95 ±0.04 


0.95 ±0.04 


0.23 ±0.01 


Total Heavy 1=0 


88.0 ±4.7 


88.0±4.7 


13.3 ± 0.7 


QCD continuum 








3. -> (4.57) 2 
(4.57) 2 -» (11.27) 2 
(11.27) 2 -> 4M t 2 
4M t 2 -> 00 


407.0 ±2.3 


407.0 ±2.3 


49.4 ±0.3 


95.3 ±0.5 


95.3 ±0.5 


27.3 ±0.1 


20.5 ±0.1 


20.5 ±0.1 


5.87 ±0.01 


w 0. 


w 0. 


W 0. 


Tote/ QCD CW. 


5£g.<5 ± 2.^ 


± 2.-4 


82.6 ±0.3 



|3C,31] 

|3gQ 



NWA ||| 
NWA || 

^exclusive (34 

BW ||,|| 

BW 



NWA 
NWA 



[33| 
[33| 



(u, d, s) 
(u, d, s, c) 
(u, d, s, c, 6) 
(it, d, s, c, 6, t) 



7036.1(76.4) 



7015.7(119.1) 



354.4(3. 



decays predominantly into pir (we assume it to be ap- 
proximately 90%); the w(1662) decays mainly into pir 
(44%) and wtt (56%), while the 0(1680) decays domi- 
nantly into K*K (93%) and to KK (7%). In this way, 
we deduce the leptonic widths: 



- u;(1419)->e+e- 
^w(1650)->e+e- 
^0(168O)^e+e- 



(90 ± 34) eV , 
(0.30 ± 0.03) keV , 
(0.43 ± .15) keV . 



(17) 



We use a Breit-Wigner (BW) form for evaluating their 
contributions: 



<*H{t)\ 



BW 



12tt- 



r„ R r* 



(t-M 2 ) 2 + M 2 rL ' 

where we take Mr from PDG |33] and: 



(18) 



w(1420) 
tot 

w(1650) 
tot 

0(1680) 
tot 



(174 ±59) MeV [35], 
(280 ±24) MeV ||, 
(150 ±50) MeV ||. 



(19) 



3.3. THE J/ip AND T FAMILIES 

We consider the six J/ip mesons IS*, 25, 3770, 4040, 
4160 and 4415, and the five T mesons 15, 25, 35, 45 
and 11020. We treat these mesons using NWA and the 
data on their electronic widths from PDG [p3|| . 

3.4. QCD CONTINUUM 

As advertized previously, we shall treat the region 
not accessible experimentally by using perturbative 
QCD where it is expected to work and where non- 
perturbative effects like e.g. the quark and gluon con- 
densates are negligible corrections. 

Input for light quarks 

In the case of massless quarks, the unit operator con- 
tribution to R e + e - is accurately known and reads, in 
the MS scheme: 



F 3 a 2 s + F A a 3 s ] 



(20) 



/ 



with: 

F 3 = 1.986- 0.11571/ [37] 
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F A = -6.637 - 1.200ra/ - 0.005ti/ 



-1.240 



(g/g/ 
3 (E/ Q 



(21) 



for rt/ quark flavours; Qf is the quark charge in units 
of e and a s = (a s /-7r) is the QCD running coupling 
which we use to order a 2 (see e.g. the expression in 
We also use, the value of the QCD scale for 
three flavours: 



A 3 = (375 ± 50) MeV, 



(22) 



obtained using the central value a s (Mz) = 0.119 [[53], 

To this perturbative correction, we add the light 
quark mass and non-perturbative corrections. We shall 
also consider the effects of a tachyonic gluon mass eval- 
uated in |^2| in order to take into account the effects 
of the truncation of the QCD perturbative series. 

The dimension D = 2 (m 2 and tachyonic gluon) 
corrections are ||,||,|3|: 

= -3^Q 2 m/[6 + 28a s + 
/ L 

(294.8- 12.3n/)a 2 ] + 1.05a s A 5 



e+ e 



, (23) 



where g§ a s X 2 ~ -(0.06 - 0.07) GeV 2 
The dimension D = 4 corrections are 



t z R 



(4) 



3J2Q 

f 



2tt / 11 

— 1 a, 

3 I 18 * 



,G 2 



+8tt 2 (1 — a s ) mf(tpftpf} 
32tt 2 — 



27 



6 23 



+mi-\ + — 

s 7 L a s 4 

/2063 



V 24 



-100 



(24) 



where Cs = 1.2020569... 

We use (a s G 2 ) = (0.07 ± 0.01) GeV 4 g|, and 
(m u + m d )(uu + dd) = -Iflml (/, = 92.6 MeV). 
The running quark mass in the MS scheme is defined 
as (see e.g. ]4p|]): 

mi{v) = 771* (-ft a^))" 71 ^ 1 |l 
. ft ( 7i 72 \ f s 



1 


[1 


f 7i 


\ 2 

72 \ 


/3 2 2 


( 7i 


72 


2 




Ui 


/32; 1 


0i 


Ui 


02 



{h fix 

"a Ui 



73 

ft 



(25) 



where for three flavours the quark mass anomalous di- 
mensions are (see e.g. pj|f4C|): 71 = 2, 72 = 91/12, 



73 = 24.8404, while the coefficients of the (3 function 
are: ft = -9/2, ft = -8, ft = -20.1198. We shall 
use the value of the invariant strange quark mass ]46| : 



m s = (133.3 ± 18.8) MeV 



(26) 



We take into account the SU (3) breaking of the quark 
condensates (ss) / (dd) ~ 0. 68^0 29 Q and consistently 
use f K = 1-2/tt. 

Input for heavy quarks 

In the case of heavy quarks, the perturbative spec- 
tral function can be accurately approximated by the 
Schwinger interpolating formula |4q|: 



R. 



o (3-«?) r 4 
!+e - =3^Q 2 w/ i_^[i+ -«./(«/)] 



where: 



/(«/) = 



and: 



7T (3 + U/) (it 

2u7 4~ 



2 4tt 



v f 



Am 2 
t 



(27) 



(28) 



(29) 



is the quark velocity; mt is the heavy quark (pole) 
mass of flavour /. We shall see that, in the region 
where we shall work (away from threshold and for large 
t), these two parametrizations provide a sufficiently ac- 
curate description of the spectral function. For a con- 
servative estimate, we shall consider the range of values 
spanned by the running and pole quark masses for the 
heavy quarks |j49]j : 

m c = (1.2 - 1.46) GeV, m b = (4.2 - 4.7) GeV. (30) 

Within the present accuracy of the values of the heavy 
quark masses, we find that it is not necessary to include 
the known a 2 corrections ijsof . 

The (u, d, s) flavour contributions 

We complete the contribution of the light quark chan- 
nels by adding the QCD continuum from t c — 3 to 
(4.57 GeV) 2 . We shall use the value of A 3 given in Eq. 
(|22|) . Our value of t c is higher than t c ~ 1.6 GeV 2 from 
global duality constraint in the 1=1 channel fl5llfj. We 
have compared a^ =1 (1.6 < t < 3) GeV 2 from a QCD 
continuum parametrization and from the uses of data. 
They only differ by 16xl0~ n , showing that the QCD 
continuum provides a good smearing of the data. We 
shall include this value into the systematics. One can 
also realize that the mass and non-perturbative correc- 
tions tend to cancel each others, whilst, individually, 
they are already small. 

The (it, d, s, c) flavour contributions 

We add the contribution of the charm quark about 150 
MeV above the J/f/)(4415) (~ empirical mass-splitting 



1 I thank Eduardo de Rafael for this remark. A similar value of 
the QCD continuum is obtained from the V- A channel |f52j. 
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between the radial excitations). In this region, the pre- 
vious expressions give a good description of the spec- 
tral function, because at this energy, the charm quark 
is already relativistic with a velocity larger than 0.75. 
Additional threshold effects are taking into account by 
transforming the value of A 3 into A4. 

The (u, d, s, c,b) flavour contributions 

The contribution of the &-quark is added about 250 
McV above the T(11020) resonance (« empirical mass- 
splitting between the radial excitations). Again at 
this energy, the b quark is already relativistic as its 
velocity is larger than 0.55, such that our perturba- 
tive parametrization remains a good approximation. 
We integrate until the 2M t threshold, where we take 
M t ~ 174.3(5.1) GeV as a best value given by PDG 

M. 



3.5. FINAL RESULT FOR 



-.had 



{1.0) 



From combined r-decay and e + e data 

Collecting all different contributions from Table 2, we 
deduce the final result at the end of that table, which 
we compile in Table 1 for a comparison with other pre- 
vious determinations. We obtain ^: 

al ad {l.o) = 7036.1(76.4) x 10~ n , (31) 

using the inclusive r-decay data until 3 GeV 2 from 
ALEPH/OPAL |^,|l| for the isovector channel and 
using e + e~ data below 3 GeV 2 J33|-|3(| for the isoscalar 
channel. The errors have been added quadratically. 
The effects of the heavy quark mesons have been 
treated using a narrow width approximation plus QCD 
continuum away from the quark- anti-quark thresholds. 

From the alone e + e~ data 

If we use the isovector e + e~ data below 3 GeV 2 , we 
obtain: 



had 



(l.o) = 7015.7(119.1) x 10 



-11 



(32) 



which is slightly lower than the one from inclusive r 
decay and less accurate. The difference in each region 
can be seen explicitly in Table 2 and easy to understand 
from the data given in [B0 24 . 



Different source of errors 

— The main error (80% when added quadratically) in 
our previous determinations comes from the p-meson 
region below 0.8 GeV 2 . Hopefully, improved measure- 
ments of this region are feasible in the near future. 

— The second source of errors comes from the re- 
gion around M T for the inclusive r-decay and between 
1 GeV to M T for the e + e~ data. These errors are 
about half of the one from the region below 0.8 GeV 2 
in most different determinations. They can be reduced 
by improved measurements of inclusive r-decay near 
M T (1=1) and by improving the measurements of 
the odd multi-pions and KK, KKtt, ... channels in the 
1 = channels from e + e~ data. 

— The contributions of the whole region above 3 GeV 2 

2 Hero and in the following, the errors will be added quadratically 
as usually done. 



induce much smaller errors (7% of the total). There 
is a quite good consensus between different determina- 
tions in this energy region. 

Improvements from r-decay data 

Our results for the isovector contributions, show that 
the precise inclusive r-decay data have significantly im- 
proved by almost a factor two the accuracy of the de- 
termination compared with the one from e + e~ data. 
However, they are amended by the inaccuracy near 
M T . Indeed, one can gain a bit in precision by us- 
ing e + e~ data in the region between 2.1 and 3 GeV 2 
instead of r-decay data, but we are aware of the dis- 
agreement among different e + e~ data in this region. 
In this case of figure, one would get from Table 2: 



.had 



(l.o) = 7010.0(72.2) x 10 



-11 



(33) 



Final result and additional systematic errors 

We consider as a final result the average of the three 
results given in Eqs. (|3l|), (|32] ) and (p3|). We consider 
the accuracy given in Eq. (|33|), which avoids the weak- 
ness of r-decay data near M T . We add to this error the 
following systematics: 

— We take into account, an eventual deviation from 
the CVC assumption (CVC has been also tested 
in j53|), and some other systematics in manipulat- 
ing the isovector data by adding a systematic error of 
15.5xl0~ n coming from the largest difference of the 
central value of Eqs. ( ^l|) and ((3^) from the average 
value. 

— We also add another 16xl0 -11 systematics due to 
the difference of the continuum threshold 3 GeV 2 value 
for three flavours used here and the one from QCD 
global duality in []5l| (see section 3.4). 

Adding the different errors quadratically, we deduce 
our final estimate: 



al ad {l.o) = 7020.6(75.6) x 10 



-11 



(34) 



It is amusing to notice that the central value of our 
different results almost coincide with our 25 years old 
value HQ and with the mean central value given in 
Table 1. 

3.6. COMPARISONS WITH PREVIOUS DE- 
TERMINATIONS OF a* ad (l.o) 
QCD continuum 

— Our treatment of the QCD continuum in the heavy 
quark channels gives similar results than more involved 
parametrization p2[ including non-perturbative gluon 
condensate and non-relativistic contributions. This is 
due to the fact that in the region where we work, 
the heavy quark-antiquark pairs are already relativis- 
tic with a quark velocity larger than 0.55. 

— Comparing existing estimates, one can see that 
there is a consensus on the size of the QCD continuum 
effect. The departures from different determinations 
come mainly from the treatment of the low-energy re- 
gion below 2 GeV. 



3 1 thank Paco Yndurain for mentioning this point, for checking 
some of the entries in Table 2 and for several communications. 
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Isoscalars 

For the isoscalar resonances, one can realize by com- 
paring the results with existing estimates that the 
narrow width approximation (NWA) gives enough ac- 
curate results comparable in magnitude and accu- 
racy with more sophisticated parametrization includ- 
ing width effects [^lj . 

Isovector channel and some global comparisons 

In the following comparison, we shall use numbers in 
units of 10 _u . 

— Comparison with JjljJI^/ 

It has been argued in |p_9|| that the contribution from 
this region can be improved by using constraint im- 
posed by 7T7T scattering data in the spacelike region 
on the pion form factor. The result a^ ad (l.o)(t < 0.8 
GeV 2 ) = 4848(31) quoted by |l| is based on Wat- 
son theorem (extended to 0.8 GeV 2 ) for the 7T7T phase 
shift S{(t). The central value differs with our value 
4795(61) obtained from the data in the timelike region 
by 53 MeV, but agrees with the result 4794(14) given 
in their Eq. (3.23) obtained from a fit of 8l(t) from the 
data. Another discrepancy of 39(104) comes from the 
region between 0.8 and 2 GeV 2 , but is not significant 
due to the large errors. It is argued in p3] that error 
from high-energy region can be reduced . 

— Comparison with J2fy 

Our result agrees within the errors with the final re- 
sult 7024(153) of H from their Table 3a 0, with the 
updated result 6989(111) quoted in [§J, and with some 
of the results given in Table 2. 

— Comparison with J2^^2^J 

For a proper comparison with these different de- 
terminations, we shall only consider the two regions 
separated by the scale \ft = M T , and compare our re- 
sult a* ad (l.o)(Vi < Mr) = 6426(76) obtained from t- 
decay data. The central value of our result is relatively 
higher than the value 6365(74) in || and 6343(60) in 
[ p6]| . Within the errors, the 1st result is consistent with 
ours while the second shows one a discrepancy. For 
comparing with the result of pi] , we subtract from 
their total result, the heavy quarks and QCD contin- 
uum contributions. In this way, we deduce the contri- 
bution 6405(94), much more closer to our result than 
the two formers, where one should also notice that our 
final results from the alone e + e~ data agree. 

The exact origin of the discrepancy of the central 
values is not easy to detect due to the complexity of 
the analysis in this region, though it is reassuring that 
results based on maximal data inputs J24|,|25|] are con- 
sistent with ours within the errors. 

We do not worry by the discrepancy with |2(J, as 
it uses more theoretical inputs than data ("determi- 
nation with minimal data input" as quoted from [ p^ ) 
for minimizing the errors , through the introduction 
(and subtraction) of some arbitrary polynomials in t 
(not coming from QCD first principles) for approach- 



ing the t — region 0. On the other hand, the use 
of "local duality" in a given interval is stronger than 
the standard applications of QCD spectral sum rules 
a la SVZ Q (for a review see e.g. 0]), where sys- 
tematic uncertainties can be more severe. We do not 
expect that the result obtained in this way is much bet- 
ter (though apparently more accurate) than the ones 
previously obtained from pure data analysis. 

One should also notice that, though one has been 
able to formulate the r-decay width using QCD and 
the Operator Product Expansion (OPE) including the 
QCD vacuum condensates §|,|| a la SVZ @, which 
has provided an accurate measurement of the QCD 
coupling a s , the extension of such a program to 

the analysis of the lepton anomaly cannot be done in a 
straightforward way due the existence of a pole at t = 
0, which needs a nonperturbative QCD control of the 
contribution of the small circle at the origin when using 
the Cauchy contour for evaluating the muon anomaly 
in the complex-t plane. Such a control is not feasible 
at present from QCD first principles. 

Final errors and comments on the uses of new 
data 

We can conclude from the previous determinations in 
Table 1 and the present work that, with present data 
and with less theoretical inputs, a realistic error for 
extracting a^ ad (l.o) cannot be less than 70 x 10~ n . 
After the completion of this work, we have checked ||(| 
by doing similar analyses (hadronic contributions to a 
and to the muonium hyperfine splitting) that, using 
the new Novosibirsk |5lj and BES |58| data, one does 
not improve noticeably the errors compared to the one 
obtained from the correlated averaged data and QCD 
continuum parametrizations used in this paper. 

4. THEORY OF THE MUON ANOMALY 

4.1. SUM OF THE THEORETICAL CON- 
TRIBUTIONS 

— QED and SM electroweak contributions 
Using obvious notations, they read PJ59|,|5||: 



QED 



EW 



116 584 705.7(2.9) x 10 
151(4) x 10 



-li 



-li 



(35) 



with QED up to 8th order and EW including two-loop 
corrections pof . 

— Higher order hadronic contributions 
For a^ ad (/i.o), we have used the contribution of the 
high-order vacuum polarizations (h.o)v.p recently es- 
timated in |3l| , where we have checked that the central 
value coincide with the original estimate in [|| which 
has the same result for a^ ad (i.o) though the error is 
larger in |6| due to older data. This feature is reassur- 
ing for a self-consistency check. As the data used in 
give smaller value of a^ ad (l.o) by 1.00125 than 



4 I thank William Marciano for mentioning this correct result 
and for several communications. 



5 This procedure together with the analytic continuation method 
has been already criticized in | M0l . Analogous comments apply 
to the estimate of a(Mz)- 
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the present work, we have rescaled the (h.o)v.p contri- 
bution by this small factor, and we deduce: 



had 



(h.o) v . P = -101.2(6.1) 



(36) 



For the light by light scattering hadronic contributions, 
we use the average of the two estimates in [ [52"| which 
agree each others. We have corrected the sign of the 
pseudoscalar and axial-vector contributions according 
to the result in [Q, which has been confirmed in |64|] 
fl. Then, we obtain: 



had 



(h.o) LL = 70.8(21.0) x 10 



-li 



(37) 



We have not included into the sum the contribution 
due to the quark loop, suspecting that adding this ef- 
fect together with the pion loop and meson pole contri- 
butions can lead to a double counting (quark-hadron 
duality arguments []). Instead, we have considered the 
quark loop effect incorporating vector mesons in the 
photon legs, which is about 15 x 10 -11 (average of the 
two determinations in |62|) as another source of un- 
cerainties of the estimate. Then, we deduce: 



had 



(h.o) 



„, h a d 



(h.o) v . P + a™(h.o) 
-30.4(21.9) x HT 11 , 



LL 



(38) 



showing that the higher order hadronic effects tend to 
cancel each others. 

— Total theoretical predictions in the SM 

Adding the results in Eqs. ( |35| ) and (^8|), one obtains 

the SM theoretical contributions: 

al M = 116 584 826.3(22.5) x 10" 11 + a£ ad (Lo) 

= 116 591 846.9(78.9) x 10" 11 , (39) 

where a^ ad (Z.o) is the lowest order hadronic contribu- 
tions evaluated in this paper [see Eq. (|34|)] . 

4.2. THEORY VERSUS EXPERIMENT 

In confronting the theoretical estimate with the ex- 
perimental value in Eq. (Q), we shall use, for reasons 
explained previously, the result in Eq. ( |34| ) with a 
moderate accuracy. We deduce: 



new _ ~exp 
a H — a n 



,SM 



176(170) x 10 



-li 



(40) 



The mean central value from our selected determina- 
tions given in Tabic 1 would give: 



,SM\ 



189(175) x 10 



-li 



(41) 



after adding to the "statistical" error in Table 1, the 
one from our determination in Eq. (|34|). Both values 
indicate an eventual window of about only a for new 
physics beyond the standard model. A more definite 
claim is waiting for a more precise et M measurement. 



°We are also aware of the fact that the two results are model- 
dependent and further tests on their validity are needed. How- 
ever, one can notice that this contribution is of the same size as 
the error of a l ^ ad (l.o). 

7 The agreement o f tin e result in_Eq. ( [j?] ) with the one 62(3) X 
10 -11 obtained in |18| (see also fej), using the quark constituent 
model (and no hadrons at all!) seems to support such arguments. 



and (El]) into 



(42) 



We translate the results in Eqs. ( pj| ) 

the conservative range: 

- 42 < a" elu x 10 11 < 413 (90% CL) 

We shall discuss some numerical implications of this 
result by giving conservative lower bounds on the scales 
of new physics discussed recently in the literature |^|, 
66|-|9). 



4.3. LOWER BOUNDS ON NEW PHYSICS 
Supersymmetry 

Using the expression given in Q and for tan/3 e 
(</>2)/(</>i) (ratio of Higgs expectation values) > 4: 

, , , 2 



-,SUSY 



a 



■ tan (3 



(43) 



one obtains a conservative lower bound on the degen- 
erate sparticle mass: 



m > 113 GeV (90% CL), 



(44) 



which is comparable with the present experimental 
lower bound of 100 GeV. 

Radiative mass, composite and extended Zee 
models 

In these models where the typical contribution to 
is about (m 2 /M 2 )0(l) one can derive thc 

lower bound on the scale of the models: 



M > 1.7 TeV . 



(45) 



This constraint can also be translated into a lower 
bound on the mass of a singlet scalar in the Zee Model 
with a coupling 1/967T 2 |p8fl. One can deduce: 



M S2 > 55 GeV 



(46) 



Leptoquarks scenarios 

Leptoquark effects on have been discussed recently 
in fi9"[ . For an electromagnetic coupling at the lepton- 
leptoquark-quark vertex, one can deduce the approxi- 
mate expression: as: 



4-0 



giving: 

M lq > 1.1 TeV 



(47) 



(48) 



This bound is much larger than the present lower 
bounds of about (200 ~ 300) GeV from direct search 
experiments at HERA and Tevatron. 
We expect that these different bounds will be improved 
in the near future both from accurate measurements 
of and of e + e~ data necessary for reducing the the- 
oretical errors in the determinations of the hadronic 
contributions, being the major source of the theoreti- 
cal uncertainties. 

5. THEORY OF THE TAU ANOMALY 

Here, we update our old and first work [Q on the t 
lepton anomaly. We shall use M T = 1.77703 GeV from 
PDG 11, 
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Lowest order hadronic contribution 

The estimate from e + e~ + r-decay is given in Table 2. 
Taking the similar average as in the case of the muon 
(section 3.5), one obtains: 



a^ ad (l.o) = 353.6(4.0) x 10" 



(49) 



One can see that the estimate 280(20) x 10~ 8 in [pf, 
using present values of the QCD condensates |^5|, 4C | , 
is much lower than this result. This appears to be a 
general feature of this method based on minimal data 
inputs (see e.g. |2(| for a^ ad (l.o) discussed in the pre- 
vious section). A comparison with the result of [^H 
shows that in both cases of muon and tau, this result 
is systematically higher than ours by the same scaling 
factor 1.0004. Previous results 370(40) |§ and 360(30) 
[fnj in units of 10~ 8 are comparable with ours but in- 
accurate. 



Higher order hadronic contributions 

The higher order contributions due to vacuum polar- 
ization a'^ ad (h.o)v.p, can be obtained from the result 
of 0. We rescale it by the factor 1/1.0083 like we 
have done for the muon case as explained in previous 
section 9. For the light by light scattering contribu- 
tion, we use the result of |6^] for the muon and we 
rescale it by the mass squared ratio (m T /m p ) 2 , which 
is expected to be a good approximation from the semi- 
analytical expression given in 0. Then, we obtain: 



a'; ad (h.o) v .p 
a h T ad (h.o) L . L 



7.6(0.2) x 10" 

2 



20.0(5.8) x 10" 



Therefore, we deduce: 
a h T ad {h.o) = 27.6(5.8) x 10" 
and: 



(50) 



(51) 



had n had i 



(l.o) + a h T ad (h.o) = 381.2(7.0) x 10" 8 . (52) 



Electroweak contributions 

The lowest order electroweak SM contribution is well- 
known: 



EW 



(l-o) = r 



5 G T ml 



1 + -(1 - 4sm 
5 



2 V) 2 



(53) 



3 8tt 2 V2 

to order (m r /M) 2 (M being the W or Higgs mass); 
We assume G T = G M = 1.166 39(1) x 10~ 5 is the 
Fermi coupling; sin 2 9 W = 0.224, where 9w is the weak 
mixing angle. We add the two-loop contribution (2- 
loop/l.o~ — 65a/7r) |30|, which induces a 15% reduc- 
tion of the one loop result. We obtain: 



EW 



46.9(1.2) x 10" 



(54) 



QED contributions 

— Generalities 

These contributions have been first evaluated in Q| to 



order a 3 . Here, we revise and improve these evalua- 
tions. In so doing, we use the relation ^,^3j: 



mi (=)■ 



1,1' =e,fi 

[a 6 (l) + a 6 



(55) 



where I indicates the internal fermion loop appearing in 
the photon propagator, (I) and a^il) can be deduced 
from the result of the muon anomaly while a§(ll') is a 
new contribution involving one electron and one muon 
loop insertion in one photon propagator. At the same 
order of truncation of PT series, one has [[73] , j59| : 

2 



a. 



1.181 241 456 (-) 

1-8 



- 0.328 478 965 

3 



+ 



= 1 159 65.2 x 10" 
where we have used the measurement: 
of 1 = 137.036 003 7(33) , 



(56) 



(57) 



from the quantum Hall effect [|74j . 
— QED at fourth order 

For evaluating this contribution, we use the dispersive 
representation: 



1 f°° dt / 4m 2 

o*(0 = | f (* + 2m 2 )^/l - R T (t) , (58) 

where K T (t) is the kernel function defined in Eq. (||). 
Then, we obtain: 

a 4 (e) = 2.024 29 and a 4 (£t) = 0.361 66 . (59) 

The former result is well approximated by the known 
analytic approximate relation to order (m e /m r ) 2 given 
e.g. in H , while the second result needs the inclusion of 
the not yet available (m p /m T ) 3 term Adding these 
two contributions, we obtain: 



(a T - a e ) 4 = 2.385 95 



1287.3 x 10" 



(60) 



— QED at sixth order 

Because of the accurate determinations of the hadronic 
and weak contributions, the inclusion of the sixth or- 
der contributions becomes necessary as they contribute 
with the same strength. The contributions of dia- 
grams with vacuum polarizations and ladders, can be 
obtained from the analytic result given for the muon 



a 6 (l) = hog 2 ^- 1.113 90 log ^ 
9 mi mi 

4.307 66 + ... , 
which leads to: 



(61) 



a 6 (e) = 10.000 2 , a 6 (^) = 2.934 . (62) 

8 These results agree with the last paper in |nj confirming that 
the results in their two former papers are wrong. 
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a§(n) differs from the value —.122 given in Jt3] where 
the origin of the negative value is not understandable 
from Eq. (|6l| ) but should come from a wrong term used 
there. We add the new class of contributions specific 
for a T given in 0: 



a 6 (e/i) = 



12 m T m 7 
2< o log log — 

9 m„ m. 



25, 

— log Z- 

54 m P m,, 




2.753 16 . 



(63) 



We parametrize the light by light scattering contribu- 
tion as: 



27T^ TYl 

a 6 (l) L . L ~— log — 1-14.13 , 

3 TO/ 



(64) 



by combining the known coefficient of the log-term ]7q ] 
and the numerical value of the total contribution for 
the muon ffql. This leads to: 



« 6 (e)L.L = 39.521 7 



o 6 (m)l.l = 4.441 2 , (65) 



This contribution is so far the most important at sixth 
order. Adding the different contributions from Eqs 
(EHl) to (p5|), one obtains to sixth order: 



(a T - a e ) 6 = 59.650 (-J = 74.758 



x 10" 



(66) 



This effect is about the same strength as the weak 
interaction effect and bigger than the higher order 
hadronic contributions. 

— Total QED contribution up to sixth order 
Adding the previous QED contributions, we deduce 

a QED _ n7 327.1(1.2) x 10~ 8 , (67) 
where the error comes from the measurement of a. 
Final result 

Summing up the previous different theoretical contri- 
butions, we deduce in the standard model: 



Q SM = a had + a EW +a QED 



117 755.2(7.2) x 10" 



(68) 



We consider this result as an improvement of our old 
and other existing results. In |^6|, only an average 
of different existing hadronic contributions have been 
added to the QED and electroweak contributions. In 
J7|]7l|], the hadronic contributions are inaccurate, while 
the value of the sixth order QED contribution is incor- 
rect. This value in Eq. ( |68| ) can be compared with the 
present (inaccurate) experimental value ff% '. 

a e T xp = 0.004 ±0.027 ±0.023 , (69) 

which, we wish, will be improved in the near future. 

6. SUMMARY 

— We have re-evaluated the lowest order hadronic con- 
tribution to the muon and tau anomalies using the pre- 
cise T-decay data below M T and averaged e + e~ data. 



Ours results are given in Table 2 and section 7, and 
our best choice in Eq. ([m|). Though the approach is 
not conceptually new, the present situation of a M (the- 
ory versus experiment) justifies a new independent re- 
evaluation of the anomaly using minimal theoretical 
inputs. 

— We have extensively discussed the different sources 
of the errors in the analysis, which are dominated by 
the region below M T . However, the relative weight 
of this region decreases for a T (see also detailed discus- 
sions in ) where a more precise theoretical prediction 
can then be provided. 

— We have compared our analysis summarized in Ta- 
ble 2 and our final result in Eq. (pj) with some of the 
existing estimates given in Table 1 PL 

— We have shortly discussed implications of our result 
to some models beyond the standard model, namely 
supersymmetry with large tan/3, radiative muon mass, 
composite, some extended Zee and leptoquark models. 
The lower bounds on the mass scale of these models 
from the muon anomaly is comparable or in some cases 
stronger than existing experimental lower bounds. 

— We have completed our work by updating our old 
estimate [Q] of the different theoretical contributions 
(QED up to sixth order, higher order hadronic and 
electroweak) to the r lepton anomaly. Our result is 
given in Eq. (|68|). A precise future measurement of 
the r-lepton anomaly is welcome as it will permit to 
probe the QED series at shorter distance, and where, 
the relative weights of different interaction effects are 
very different from the case of the muon. 
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